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We examine string (vortex) formation at a quench for a weakly-coupled global U(l) theory when 
the excitation spectrum is non-relativistic. It is so similar to vortex production in the corresponding 
relativistic plasma as to reinforce arguments for the similarity of vortex production in the early 
universe and in low-temperature many-body physics. 
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_— , . For many years it has been argued by Kibble p| and others M that the large-scale structure of the universe 
I ' can be attributed to cosmic strings formed during a symmetry-breaking transition at the Grand Unification scale. 
Unfortunately, given the unlikely event of observing a cosmic string (vortex) directly, chains of inference are sufficiently 
' long that it is difficult, if not impossible, to make the case compelling. 

However, the formation of topological defects like vortices during symmetry-breaking phase transitions is not unique 
to the early universe but generic to many physical systems. In particular, recent experiments on the production 
of vortices in superfluid 4 He |jj and 3 He || have excited considerable interest. The suggestion that these, and 
other experiments || , may provide an insight on the early universe have been made by the experimental grou ds 
concerned. In this they have been championed by some theoretical astroparticle physicists, most notably Zurek (| 
, whose most recent review article M was, indeed, titled Cosmological Experiments in Superfluids and Superconductors. 
• Nonetheless, despite the hard thinking that has already taken place, there is no doubt that a fuller understanding of 
the nonequilibrium dynamics of vortex production is required if comparisons are to be more than superficial analogies. 

What encourages us in the hope that vortex production in the early universe and the laboratory has close parallels 
is that, in the first instance, the most plausible production mechanism refers to neither. We recapitulate it now for the 
simplest theory permitting vortices, that of a complex scalar field </>(x, t). The complex order parameter of the theory 
' is {4>) = ne la and the theory possesses a global 0(2) (or U{1)) symmetry that we take to be spontaneously broken at 
a phase transition. The field (f> could be cither a complex non-relativistic order field appropriate to a superfluid or a 
relativistic field in the early universe. The transition is continuous in both cases but, as we shall see, the details of 
D , the transition order are largely irrelevant to our conclusions, as long as it is not strongly first-order 
-^3 • Initially, we take the system to be in the symmetry-unbroken (disordered) phase. We have no reason to choose 
any particular initial field configuration, beyond the requirement that the field is distributed about <f> = with zero 
mean r\ — \{$)\ — 0. The simplest assumption is that, beginning at some time t — to, the 0(2) symmetry of the 
ground-state (vacuum) is broken by a quench, a rapid change in the environment inducing an explicit time-dependence 
in the field parameters. Once the quench is completed the 0-field potential V(cf>) — — a\<j)\ 2 + b\(j>\ 4 is taken to have 
the symmetry-broken form a > 0, b > of the familiar 'wine-bottle' bottom. The ground-state manifold (the circle 
S 1 , labelled by the phase a of ((f))) is infinitely connected and the theory possesses global strings or vortices, labelled 
by a winding number n £ Z. Specific straight-string solutions to the classical field equations are well-known, both 
for the non-relativistic || and the relativistic cases The only property of these solutions that we need is that, as 
tubes of false vacuum, their thickness is the Compton wavelength of the massive (Higgs) excitations of the theory. 

These strings cost considerable energy to produce. That they should appear at all follows from a general argument, 
due to Kibble (IJ, which goes as follows. During the transition, the complex scalar field begins to fall from the false 
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1 ln that case the mechanism for the transition, bubble nucleation, would lead to very different consequences from those 
outlined below. 
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ground-state into the true ground-state, choosing a point on the ground-state manifold at each point in space, subject 
to the constraint that it is continuous and single-valued. For continuous transitions, for which this collapse to the 
true ground-state occurs by spinodal decomposition or phase separation, the resulting field configuration is expected 
to be one of domains within each of which the scalar field has relaxed to a constant ground-state value (i.e. constant 
field magnitude and phase). 

If this is so, then the requirements of continuity and single valuedness will sometimes force the field to remain in 
the false ground-state between some of the domains. For example, the phase of the field may change by an integer 
multiple of 2tt on going round a loop in space. This requires at least one zero of the field within the loop, each of 
which has topological stability and characterises a vortex passing through the loop. The density of strings is then 
closely linked to the number of effective domains and the evolution of this density is, correspondingly, linked to the 
nature of the domain growth. When the phase transition is complete and there is no longer sufficient thermal energy 
available for the field to fluctuate into the false ground-state, the topological defects are frozen into the field. Simple 
counting arguments suggest defect densities at this time of the order of one string passing through each correlation 
area. From then on, the defect density alters almost entirely by interactions of defects amongst themselves, rather 
than by fluctuations in the fields |J, but such calculations go beyond the scope of this paper. 



In an earlier paper [|10[, henceforth known as /, (but see also jyj), two of us (A.G and R.R) showed how defects were 
produced in a very simple model of a relativistic 0(2) scalar field theory in which the long wavelength fluctuations 
that drove the transition were Gaussian. The fluctuations were taken to be set in train by the implementation of a 
quench from a high temperature, above the transition, to a low temperature. Despite its oversimplifications, its main 
interest is that it provides a concrete example in which the scenario outlined above is, indeed, true and quantifiable. 

With low temperature many-body theory in mind, in this paper we show, in the same 0(2) model of Gaussian 
fluctuations, how a density quench of a nonrelativistic medium also leads to vortices. A priori, we might expect 
significant differences in the vortex production. The relativistic regime considered in / was characterised by initial 
temperatures T ^> m for all particle masses m (in units ks = c = 1). On the other hand, non-relativistic media 
are characterised by T <§; to, to freeze out antiparticles, but with T ~ fj, m , the non-relativistic chemical potential. 
Surprisingly, in our simple model we find no difference once a proper identification of mass parameters has been made. 

As it stands our model, being weakly coupled (in effect, with zero-coupling for short times), is far too simple 
to mimic superfluid helium. Equally, unless we are considering extremely weak couplings of the strength invoked 
in inflationary models, our approximation will only have limited applicability to hot quantum fields. However, in 
principle, if not yet in practice, we know how to incorporate stronger interactions within our approximation scheme. 
The identity of the weak-coupling models should be reflected in similarities in the more realistic stronger coupling 
counterparts. As a result, our main conclusion about the similarity of vortex production in the two energy extremes 
may have greater generality, even if the detailed results given here are too limited. 

The organisation of this paper is as follows. In / we considered the production of both 0(3) monopoles and 0(2) 
vortices. Since monopoles are somewhat easier to manage we concentrated on them at the expense of vortices, for 
which we quoted only such properties as were needed. In the next section we repair this omission by showing the 
way vortex densities, and density correlations, can be inferred from field configuration probabilities through field 
correlation functions. As in /, this summary is largely a recapitulation of some early work of Halperin ]T^| , modified 
to our purposes. 

One of the difficulties of making straightforward comparisons between early universe physics and many-body physics 
is the difference in formalism, a reflection of the freezing out of antiparticles in the non-relativistic regime. We evade 
this problem by showing how to incorporate both regimes in a common relativistic expression, in which the non- 
relativistic limit is obtained by tuning the chemical potential. This permits the relevant field correlation functions to 
be calculated in the same approximations as for the relativistic case. 

In the final sections we use these field correlation functions as input for the vortex density functions as shown 
in section 2. We conclude with a discussion of some of the results, and consider their implications for numerical 
simulations. In particular, the Kibble mechanism as presented above says nothing about the fraction of strings that is 
produced in small loops against the fraction of 'infinite' string (i.e. string that does not self-intersect). Infinite string 
is necessary for the creation of large-scale structure in the early universe. We shall suggest that the fraction of string 
in small loops is substantial. 



Consider an ensemble of systems evolving from a given disordered state or, more realistically, from one of a set 
of disordered states whose relative probabilities are known, to an ordered state as indicated above, producing 0(2) 
vortices. If the phase change begins at time to then, for t > to, it is possible in principle to calculate the probability 




II. VORTEX DISTRIBUTIONS FROM FIELD DISTRIBUTIONS 
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Pt[&] that the complex field 0(x, t) takes the value ^(x) at time t. Throughout, it will be convenient to decompose $ 
into real and imaginary parts as $ = ^75(^1 + ^2) (and cf> accordingly). This is because we wish to track the field as 
it falls from the unstable ground-state hump at the centre of the potential to the ground-state manifold in Cartesian 
field space. It is trivial to reconvert the real fields into complex fields (and similarly for conjugate fields, as we need 
them). 

The calculation of pt[$] will be performed later in our simple model. For the moment, consider it given. The 
question is, how can we infer the string densities and the density correlations from p 4 [$]? That we can calculate them 
at all is a consequence of the fact, noted earlier, that the string core is a line of field zeroes. This is equally true 
for both relativistic and non-relativistic 0{2) theories ^|. The zeroes of $ a (a— 1,2) which define the vortex positions 
form cither closed loops or the 'infinite' string mentioned in the introduction. Following Halperin |l2| we define the 
topological line density p(v) by 

p(r) = Yl fds^5 3 [r-K n (s)]. (1) 

n 

In (2.1) ds is the incremental length along the line of zeroes R„(s) (n—1,2,.. .) and ^Sa i s a un it vector pointing in 
the direction which corresponds to positive winding number. Only winding numbers n — ±1 are considered. Higher 
winding numbers are understood as describing multiple zeroes. If dAj is an incremental two-dimensional surface 
containing the point r, whose normal is in the jth direction, then Pj(r) is the net density of strings (i.e the density of 
strings minus the density of antistrings on dAj). 

Ensemble averaging (F[$])t at time t is understood as averaging over the field probabilities pt[&] as 

{F[$\) t = J X>$pt[*] F[$\. (2) 

We stress that, in general, this ensemble averaging is not thermal averaging since, out of equilibrium, we have no 
Boltzmann distribution. We shall only consider situations in which 

(ft(r))t = 0, (3) 

That is, we assume equal likelihood of a string or an antistring passing through an infinitesimal area. However, the 
line density correlation functions 

C ij {v;t) = {p i {r)p j {0)) t (4) 

will be non-zero, and give information on the persistence length of strings. 
It follows that, in terms of the zeroes of Q(r), Pi(r) can be written as 

Pi (r) = < 5 2 [$(r)]e ufc ^$ 1 (r)5 fc $ 2 (r), (5) 

where <5 2 [<I>(r)] = <5[$i(r)]($[<I>2(r)]. The coefficient of the 5-function in (||) is the Jacobian of the transformation from 
line zeroes to field zeroes. It permits us to define a further line density that we shall also find useful, the total line 
density p(r) 

Pi(r) = 5 2 [$(r)]|e iifc ^$x(r)9 fe $ 2 (r)|. (6) 

Unlike the case for Pi(r) 

n(t) = (pi(r)) t > (7) 

and measures the total string density in the direction i, without regard to string orientation. The isotropy of the 
initial state guarantees that n(t) is independent of the direction i. We note that the Jacobian factor multiplying the 
field (^-functions in (||) and guarantees that random field zeroes with no vorticity will not be counted. 
In general, the best that we can do is write (pi(r)) t as 



2 However, this is not necessarily the case for the more sophisticated vortices of 3 He, for example, for which different methods 
would be required. 
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Va (\e ijk d j $ 1 (r)d k $ 2 (r)\e i f dxa «**) t (8) 

Our simple model assumes that pt[ ( ^ > ( r )] i s Gaussian. Details as to why this could be will be given later, but if it is 
so then n and Cy are easily calculable. Given that 



suppose that 



that 



(3«(r)) t = 0, 
(* (r)a,$ 6 (r)) t = 0, 



($„(r)* 6 (r , )) t = W o6 (|r-r , |;t) 
= 8 ab W(\r-t'\;t), 

and all other connected correlation functions are zero. 

Then all ensemble averages are given in terms of W[r\ t), where r = |r|. In particular, (pi(r)) t separates as 

(pi(r)) t = (S 2 Mr)]) t (levkdfrtfdkihWDt 

or, equivalently, from (||) 

It follows fl2| |, on first performing the a integration in the second factor, that 



W"{0;t) 



W(0;t) 

where primes on W denote differentiation with respect to r. Thus 
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where ^ is the length at time t that sets the scale in W(r; t). 
On decomposing the density-density correlation functions as 



Cij{T;t) = A(r-t)6 tJ + B(r;t) 



(9) 
(10) 

(11) 
(12) 
(13) 

(14) 
(15) 



(16) 



then, in the same Gaussian approximation, A and B can also be calculated in terms of W and its derivatives, with 
rather more difficulty, as 



A = 



W 



27r 2 rA 2 



(w"a + {w'yw) 



and 



in which 



IY2 



B = 



m 

27r 2 r 2 A 



A(r;t) = W(0;t) 2 - W{r;tf 



(17) 



(18) 



(19) 



We note that both the density n(t) and the correlation functions Cy (r; t) are independent of the scale of W. The 
reader is referred to fl2l for further details. 
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To understand what A(r; t) and B(r\ t) measure, suppose x = (0, 0, z). Then Cy is diagonal, with non-zero elements 

C 33 = A, C n =C 22 = A-B (20) 

For the sake of argument, consider the idealised situation in which there is exactly one string passing through each 
area n(i) _1 in the 1 — 2 plane and that these areas form a regular lattice of cell-length £. Normalise dj, A, B with 
respect to n 2 , as 

Cy = a= — , b= — , (21) 

rt z n z n z 

whence cn = a — 6, etc.. The situation in which strings passing through adjacent faces of this lattice are oppositely 
oriented (i.e. string-antistring) then essentially corresponds to cn(£;i) = C22(£;£) = —1, a case of maximum anti- 
correlation. On the other hand, if there is equal likelihood of the next face containing a string or an antistring then 
c n(£; *) = c 22(£; t) = 0. As for C33/71 = C33, this takes the value czz{£,;t) — 1 if the string is guaranteed to continue in 
the same 3-direction through the next lattice cell, and takes the value C33(£; t) = —1 if it changes direction at intervals 
£. B itself can be isolated by observing that, if x = ^(r, r, 0), then C12 = \B. Thus a value 6(£;t) = 1 is also a 
guarantee that the string changes direction every £. As we shall see, the situation is more complicated, even in our 
simple model. Nonetheless, we shall interpret strong anticorrelation in the diagonal cn (and positive b) on scales £ as 
signifying a persistence length (the typical length along the string before it has completely changed direction) of £. 

Nothing that we have said so far discriminates between vortices in a relativistic or a non-relativistic medium. This 
distinction appears in the definition of to which we now turn. 



III. FIELD DYNAMICS WITH A CHEMICAL POTENTIAL 



In order to compare vortex formation in relativistic and non-relativistic media it is convenient to develop a single 
formalism that encompasses both. For simplicity we assume in each case that the system is initially in equilibrium in 
the disordered phase and that the transition to the ordered phase occurs as the result of a 'quench', a rapid change in 
the environment of the system. The interpolation between relativistic and non-relativistic regimes is then effected by 
introducing the chemical potential fi, coupled to the conserved charge Q arising from the 0(2) symmetry. Provided 
(i is small in comparison to a particle mass the introduction of such a potential will have little effect on a phase 
transition for the relativistic theory initiated by quenching the system from a high temperature To to, effectively, zero 
temperature. However, on increasing (i prior to the transition it becomes more costly to produce antiparticles and, 
if the initial temperature Tq is decreased to a value much less than /i, antiparticles are frozen out. The system is 
then one of non-relativistic particles at a temperature much less than the particle rest masses. In this nonrelativistic 
regime the transition is induced by a quench in /i itself or equivalently, since [i determines the density, by a density 
(pressure) quench. 

The situation is summarised in Fig.^, which shows the (equilibrium) phase structure of the global 0(2) theory. 
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FIG. 1. Equilibrium phase structure of the global 0(2) theory. The arrows show the directions in which the quenches are 
performed. 



5 



The inner sector in the T — fi quadrant is the ordered phase, the outer region the disordered phase, and the 
line separating them the phase boundary. T c is the critical temperature at zero chemical potential and /i c is the 
critical chemical potential at zero temperature. The relativistic quench discussed in / corresponds to a transition at 
zero (small) fj, across the phase boundary at T c . The non-relativistic quench with which we shall be comparing it 
corresponds to a transition at small T across the phase boundary near fj, c . 
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FIG. 2. The closed timepath contour C+ © C_ . 

From this viewpoint the chemical potential is seen as determining the initial conditions for the subsequent dynamics, 
for which we adopt the closed time path method (Schwinger-Keldysh formalism) generalising the analysis 

begun in I for the relativistic theory. As a starting point suppose that, at the initial time £ ; we are m a disordered state 
with ((f)) = 0. Our ignorance is parametrised by the probability distribution p* [$] that, at time to, (f>(to,x.) = $(x). 
For the moment we take it as given. Whether we are in a relativistic or non-relativistic regime is largely encoded in 
Pt [$]. The subsequent, essentially generic, non-equilibrium field evolution is driven by a change in the environment. 
Specifically, for t > to the action for the field is taken to be 
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S[4>] = I d* X ( -d^ad^a - ^-m 2 (t)4>l ' - A mi€f 



(22) 

where m(t), X(t) describe the evolution of the parameters of the theory under external influences, to which the field 
responds. As with $, it is convenient to decompose <j> in terms of two massive real scalar fields <p a , a = 1,2 as 
= (0i + i(t>2)/"j2, in terms of which S[<j>] shows a global 0(2) invariance, broken by the mass term if m 2 (t) is 
negative. 

The change of phase that begins at time t will, by the mechanism indicated earlier, lead to the appearence of 
vortices. We saw in the previous section that the vortex distributions at later times tf > t can be read off from the 
probability [<!>/] that the measurement of 4> will give the value The evolution of pt[&] from to to tf is most 
simply written as a closed time-path integral in which the field 4> is integrated along the closed path C+ © C_ of Fig. 2, 
where cp = </>+ on C+ and <p — on C— ■ 

If T>4>± = n<t=i £^±,1 an d spatial labels are suppressed then 



pt,[9f] = I VQptM [ Vc^+V^ S[<j>+{t f ) - $/] exp(i(s[0 + ] - S[ct>-} 

J j0 ± (t o )=$ L V 



(23) 



where 8[<j)+(t) ~ $/] is a delta functional, imposing the constraint <^+(x, t) = $/(x) for each x. This is no more than 
the statement that, for a given initial state, the probability amplitude is given by the integration along C+, and its 
complex conjugate (which, when multiplied with it, gives the probability) is given by the integration back along C_. 
The ± two-field notation is misleading in that it suggests that the <fi+ (= <fi a ,+) and 4>- fields are decoupled. That 
this is not so follows immediately from the fact that 0+(t/) = (j)-(tt). 

To return to the initial conditions, we said that we could only achieve the simple analytic results of the previous 
section if pt[$] were Gaussian and, therefore, that Pt [ ( I ) ] itself be Gaussian. The simplest such distribution has $ 
Boltzmann distributed at time to at a temperature of To = Pq 1 and chemical potential fi according to a /ree-field 
Hamiltonian Ho, which we take as 



H = / d 3 x 



(24) 



where 7r a = <f> a . 

We shall motivate the use of this free-field Hamiltonian in specifying the initial field distribution later, when we 
implement our Gaussian approximations. However, since our comments on extracting the non-relativistic limit from 



G 



a relativistic field theory have greater applicability than our particular model we would like to be more general and 
extend Hq to include interactions as 



d A x 



,(V0a) 2 + ^ + jA(^) 2 



(25) 



without enforcing a Gaussian straightjacket. We stress that Gaussian does not necessarily mean free. Most simply, 
particles need to interact before they can equilibriate. More importantly, the free-field Gaussian approximation 
adopted in / and here can be extended to include interactions self-consistently in a Hartree approximation p2| . This 
will be considered elsewhere. 

Since chemical potentials are not usually relevant to relativistic bosons a few words are in order. The 0(2) invariance 
of Hq leads to a conserved Noether current with conserved charge 



d X (<p2^1 - 7T2C 



d 3 X (4>2(t>l ~ 4>2^>l)- 



(26) 



(in Minkowski space). The numerical value of Q is the number of particles minus the number of antiparticles |^. The 
thermal probability distribution p to [$] is thus taken to be 



p to [$]=($,t |e- /3o(H °-^ ) |$,to). 



(27) 



There are two ways to proceed. The first, which to us is the most natural, accepts Hq as determining the temporal 
evolution of the physical fields, and relegates fj, to the boundary conditions on the fields. From this viewpoint, Pf [$] 
can be written as the imaginary-time path integral 



pt m 



B„ffl 



2?0exp< i5"o [0] 



for a corresponding action 



So[<t>] 



d^x 



\{dAa){d^ a )~\m^l-\\( < Plf 



(28) 



(29) 



The label zero on Sq[4>] (and Hq previously) is a reminder that all are defined at time to (and not that the theory 
is free). The boundary condition B^ty] incorporates the chemical potential. In terms of the eigenstates of Q, the 
complex field <ft = (<j)i + i(f)2)/V2 and its adjoint, B^ becomes 



B^d>} : 0(to) = * = 



— (f> - e -0oM 



(30) 



where q = 1 is the 0- field eigenvalue of Q and the time-integral in (|29j) is taken in imaginary time from to to £o — i(3$. 
For (f>* , q = — 1 . In terms of the same complex fields we have 

S {<f>} = j 'dM(^0*)(d>)-mV0-A(0>) 2 ] (31) 

and, with 7r = (ni — i7T2)/v / 2, 

H a = J d 3 x[TT*n + (V0*)(V0) + m 2 0> + A(0» 2 ] 
However, for calculational purposes a decomposition in terms of 4> a , it a is usually preferable. 



(32) 



3 For a relativistic theory the natural next step would be to gauge the 0(2) or U(l) symmetry by the introduction of the 
electromagnetic field, whereby Q becomes proportional to the electric charge of the system. With non-relativistic media in 
mind, our constituents are neutral (e.g. He atoms) and we shall not do this. 
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Although we are just setting an initial condition the effect is, inevitably, to give an action So[<j>] of the form of S[<j>] 
of (^). This permits the interpretation that the action S[<j>] is valid for all times t, with the proviso that the system 
is in thermal equilibrium for t < to, during which period the mass m(t) takes the constant value m and A(t) = A, also 
constant. 

On relabelling the integration variable (j> °f (psh by $3, we now have the explicit form for p tf [$/]:- 

p tf [®f]= V$ T>(j) 3 e iSo ^ [ V^> + ^V^e i( - s ^- s ^6[(t) + (t f )-^ f ] 

J ie,!*] J <t> ± (t )=<t> 



VfaVt+Vt-. exp|iS o [0 3 ] + i(S[</>+] - 5[^_])| 5[4>+(t f ) - $/], 
where the boundary condition is now (in terms of the field combinations <$> = (<f>i + i(f>2)/V2) 

More succinctly, p^, [$/] can be written as the time ordering of a single field doublet:- 

Pt,[*f]= f V<Pe* s ^6{Mtf)~®fl 



Ba 



(33) 



(34) 



along the contour C — C+ C_ C3 of Fig. 3, extended to include a third imaginary leg, where <fi takes the 
values 0+, 4>- and 03 on C+, C_ and C3 respectively, for which 5*0 is S^+J, S[(f>-] and S 1 ^], for which last case 
m(t) = to, A(i) = A, t e C 3 . 



t - iA) 



C 4 



FIG. 3. A third imaginary leg 



As a final step in these formal manipulations we see that expression (Q) enables us to write the $-field ensemble 
averages (...)t in terms of the 0-field thermal Wightman functions. Consider the generating functional:- 



where J j<j) is a short notation for:- 



j4>= I dt [j + (t)<f>+(t) - j.<f>-(t)] + / j 3 (t)fa(t)dt, 

Jo 



-i/3 



(35) 



(36) 



omitting spatial arguments. Then introducing a a (x) where a — 1,2, we find:- 

pt,[$] = y exp|^ c [0]| expji ^ d 3 xa Q (x)[0 + (^, x) - $(x)] a | 
= J Va expj-i ^ a Q $ Q | 0, 0], 



(37) 
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where a is the source a(x, t) — a(x)8(t — tf). As with D<p, Va denotes Yii T^oi a . 

Ensemble averages are now expressible in terms of Z^. Of particular relevance, W a b(|r — r' | ; t) = ($ a (r)$b(r'))t is 
given by 

Wab{\r-r'\\t) = j 2?$ $ a (r)$ b (r') j Va expj-i J a a $ a | ^[a,0,0] 



Va 
Va 



S 2 



6a a (r)Sa b (r') 
S 2 



V$ exp 



Ua®a > Zja,0,0] 



Sa a (r)Sa b (r 



-j<S 2 [a] ^[S,0,0]i 



(38) 



On integrating by parts 



W ab {\v-T'\-t) 



Sa a {r)Sa b (r') 
(6,(r,i)&(r',t)>, 



77^ [a, 0,0] 



a=0 



(39) 



the equal-time thermal Wightman function with the given thermal boundary conditions. Because of the time evolution 
there is no time translation invariance in the double time label. 

Not surprisingly pt[$] can only be calculated explicitly in very simple circumstances, like our Gaussian approxima- 
tion, but before we do that we still have to extract the relativistic limit. Since the chemical potential is embedded in 
the equilibrium boundary conditions and, like temperature, can only be defined in equilibrium, this is essentially an 
equilibrium problem and, for the moment, we forget the dynamics. 



IV. MANIPULATING THE CHEMICAL POTENTIAL 



The expression (|34j) is valid for all /i and all To but, as it stands, is not sympathetic to the isolation of a non- 
relativistic limit. Nonetheless, the extraction of a non-relativistic regime from it is not difficult. Instead of working 
with pt [$] directly, we can work with the partition function Z M , as we saw in the previous section. The partition 
function for this equilibrium theory is (for doublet sources j a restricted to the C3-contour) Z M [0,0, 33] of (|35|), written 



Z M \fl = Z M [0,0, j] = j B V <t> expjiSbM +* j j a 4>aj, 



(40) 



where Sq[4>] is given in (E3) and we integrate only along the contour C3. As before, V(f> = VfiiVfo. On rotating to 
Euclidean time r = it, Z^[j] becomes 



Z n\j] = < Cx p| J Ja^aj) = J^V(j) eXp|-5 , S [^] - J ja^aj, 



where So, £;[</>] is the relativistic Euclidean action 



5q, 



Pa 



dr / d 6 x 



(41) 
(42) 

:,r) = e-^(x,r- i 9o) 
(43) 

where H is given in (p5|). Because of its time- independence we have dropped the t-suffix. The net charge Q (excess 
of particles over antiparticles) in the system is obtained by differentiating ZJj] — (e$ ^) as 



and the sum is taken over fields <fi = (<fii + 1^2) /v2 satisfying the boundary condition B^: 
in imaginary time. The dot now means differentiation with respect to r, cj> a = d T <fi a . 
The bracket (...) here denotes the thermal average 

(F[<j)]) =Tr{e- MH °-^ F[<j)]}. 



Pa ofi 



(44) 



3=0 
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This determines [i in terms of the chosen Q and To. In a non-relativistic environment in which antiparticles are not 
present, Q becomes the mean particle number (and hence is proportional to the density). 

Zfj\j\ suffers from the presence of the chemical potential in the boundary conditions, which makes it difficult to 
identify the phase of the system readily. This is clarified by adopting a second approach to chemical potentials, in 
which we transfer the chemical potential term [iQ to the Hamiltonian, to create an effective Hamiltonian Hq — \xQ. 
Ho — [J-Q no longer generates time translations of the physical fields <p a or, equivalently, <f>(t) and <fi*(t). However, it 
does generate time translations of the effective fields <p(t) = e %mt <\>{t) — e vt <j>(t) and <j>*(t), with q = —1. It follows 
that 4>(t) and (j)* (t) integrations, and hence <f> a integrations, are now taken over periodic configurations with boundary 
conditions Bq : 4> a (x, r) = a (x, r — /3q) in imaginary time, with no /x-dependence Details are given in the Appendix. 

The advantages of removing /i from the boundary conditions become apparent when we express of ( ffO[ ) not 
as a sum over (/)-field histories but as a sum over histories in the 0-fields. Although not usually posed this way, the 
details are well- understood. For example, see Kapusta [ p3[ , Haber & Weldon (^4|, and Bernstein et al. |25|]. On direct 
substitution, Z^[j] takes the form 



Zfilj] = J T>4> exp^-S [$',n] - J io^aj, 



(45) 



where So[<p; yt\ is now the effective action 



So[4>] lA = dr j d"x 



(46) 



where 



(47) 



and, as above, the sum is taken over periodic field configurations B , rather than the boundary condition B^. Details 
are given in Kapusta |23| and pBfl . We note that <j> a and 4> a have the same zeroes and hence are equally good for 
the calculation of vortices. It is a separate exercise to rewrite Z^[j] in terms of the conventional actions for a non- 
relativistic order field tj)(t). Details will be given elsewhere |17 but, in the interim, the basic idea is given in recent 
conference proceedings |jl|[nj by one of us (R.J.R). 

This displacement of /x from the boundary conditions to the action has enabled us to replace the classical potential 
by an effective potential in which the state of the system is more transparent. The relevant quantity is TOq, rather 
than m 2 and /i 2 separately 0. Semiclassically, when m§ > (i.e. /i 2 < to 2 ) the 0(2) symmetry is unbroken. This 
describes our disordered initial state. However, once TOq < (i.e. fj 2 > to 2 ) the free theory is unstable. This is a 
signal that the 0(2) symmetry is broken and a transition to an ordered phase has occurred. Thus, effectively, it is 
TOq, rather than to 2 that carries the time dependence, changing from positive to negative at t = to- 
Let us first consider the relativistic regime in which the initial environment is very hot, with T » m » (i. In this 
case, with \i irrelevant, the symmetry is broken by a change in to 2 . Suppose that, on the completion of the transition, 
m§ ~ to 2 takes the value Too = —M 2 < 0. If, in this relativistic case, the final temperature is very low, there are then 
no thermal effects and M is a physical parameter, determining the Higgs mass, to# = \/2M . In order to recover the 
phase boundary in Fig.l in this case it is necessary to include thermal radiative effects. As it stands too (t < 0) = too 
is not a physical parameter, but there is no loss in taking it as the effective scalar field mass at temperature To = Pq 1 ■ 
That is, in the mean-field approximation we take 



T 2 



(48) 



where To is greater than the transition temperature T c , given by T 2 — 3M 2 /A in the same approximation. It might 
be objected that this is inconsistent in that this mass is defined in terms of fluctuations at scales much larger than 
the typical domain size. Since results will turn out to be largely independent of to provided it is comparable to M 
this is not really a problem. With the parametrisation above this would not be true only if we quenched from very 
close to the transition, and we do not consider this possibility yet. 



4 We shall see later that, in our approximations, the other term depending on /i, ifx(cj)2<f>i — 4>24>i), plays no role. 
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On the other hand, for the non-relativistic theory, for which T/m <C 1, the parameter m can always be identified 
with the boson mass. We change the sign of to 2 , in the action by a change in /i from fi 2 < to 2 to /i 2 > to 2 , a density 
quench. Equivalently, in terms of the non-relativistic chemical potential 

fJ-nr = \i — m (49) 

we go from negative \i nr to positive \i nr . To preserve uniformity of notation, we take TOq (t) = m§ when t < to and 
TO oW = —M 2 < once the transition is complete. 



V. VORTEX FORMATION 



Having established the role of the chemical potential in letting us interpolate between relativistic and non-relativistic 
regimes in the initial conditions, we are now in a position to determine the effect of these initial conditions on vortex 
production in a simple model of Gaussian field fluctuations. 

We begin by recapitulating the methods adopted in / for a relativistic regime, for which /j, = 0. We have already 
assumed that the initial conditions correspond to a disordered state in equilibrium for t < to. Specifically, in I 
we adopted an initial Gaussian distribution of fields arising from the quadratic Hamiltonian Hq of (p4|). Although 
quadratic, the choice of mo > of (Es|) has thermal interactions encoded within it. We assume that interactions are 
sufficiently weak that, having established equilibrium, their effect on the Gaussian distribution is small. 

For t > to, the system is forced to change. The simplest assumption, made in I, was that, for t > to, m (i) 
and X(t) could be taken as constant, an idealised quench in which m 2 (t) changes sign to take the negative value 
m 2 (t) = —M 2 < immediately. That is, the potential at the origin has been instantaneously inverted everywhere, 
breaking the global 0(2) symmetry. Given the initial high temperature To, this can be thought of as a temperature 
quench. If A(i) = A is very weak then, for times Mt < ln(l/A), the </>-field, falling down the hill away from the 
metastable vacuum, will not yet have experienced the upturn of the potential, before the point of inflection at 
(\4>\) = y/2M 2 /3A. Thus, for these small times, X(t) can also be set to zero, and Pt[$] is Gaussian, as required. 

The onset of the phase transition at time t — to is characterised by the instabilities of long wavelength fluctuations 
permitting the growth of correlations. Although the initial value of (<f>) over any volume is zero, we anticipate that the 
resulting phase separation or spinodal decomposition will lead to domains of constant (0) phase, whose boundaries 
will trap vortices. 

We now turn to the case of non-relativistic vortex formation, cast in as similar a way as possible. The situation 
we have in mind is the following. Consider a gas of non-relativistic bosons of mass to in thermal equilibrium at 
temperature To <C to. The state is disordered. 

By changing the density we then force fi nT to change sign, from /xo < to fi al = fif > 0, as rapidly as possible. Yet 
again, in this idealised model we suppose that this change is implemented instantaneously everywhere at time t = to- 
This sets in motion a change to a superfluid phase, in which the 0(2) particle symmetry is broken by the condensate. 
Order is again established by the growth of long wavelength fluctuations in which domains form. These domains will 
trap vortices on their boundaries. 

So as to give a solvable Gaussian theory, for which we can use the results of Section 2, we maintain our Procrustean 
approach of / and ignore interactions before and after the change in chemical potential. Our initial condition is now 
equivalent to the statement that, for t < to, the particles are distributed according to the Boltzmann distribution 

n(e) = 5 , 1 \ > (50) 

where e = k 2 /2m and /xo < is the non-relativistic chemical potential \i — m. 

In the previous section we showed how the model is that of the relativistic theory of the field </> with effective action 
Sq[4>; a*] of (p6|), in which \fi n r \ ^ m - hi the notation of that section, when \i ~ m, to§ changes from 

TOq (t) = m 2 , = — 2m/io > 0, (51) 

when t < to, to 

TOq (t) = -M 2 = -2mnf < 0, (52) 

when t > to- Henceforth, we take to = 0. 

Of course, instantaneous change is physically impossible. Consider small amplitude fluctuations of <p a , at the top 
of the parabolic potential hill. Long wavelength fluctuations, for which |k| 2 < M 2 , begin to grow exponentially. If 
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their growth rate fi(fc) = \JM 2 — |k| 2 is much slower than the rate of change of the environment which is causing 
the quench, then those long wavelength modes are unable to track the quench. It will turn out that the time-scale 
at which domains appear in this instantaneous quench is td = 0(M^ 1 ). As long as the time taken to implement the 
quench is comparable to td and much less than tf — 0(M~ ln(l/X) the approximation is relevant. 
We note that, in the non-relativistic regime, f2(fc) has the same definition, but can be rewritten as 



ft 2 (k) = M - |k| 2 = 2m[n f 
= 2m((if — e(k)). 



2m 



(53) 



Thus the momentum restriction |k| < M is just e(k) < fif. 

We are now in a position to evaluate p t [$] , or rather W a b{r; t), for t > 0, and calculate the defect density accordingly. 
Details are given in the Appendix. Before we quote the result we note that the z/z(02</>i — 4>i4>\) term in So[4>; fj] 
of ( [4(f ) couples the a = 1 and a — 2 fields 4> a together and, in general, G a b{r — v';t,t') — {(j) a {r, i)</> a (r', t')) is not 
diagonal in the 0(2) labels. However, for equal times diagonal behaviour is restored as G a b(T;t,t) — 5 a bG(r;t,t) and 
the i/x(020i — 020i) term can effectively be discarded. This leaves us in a situation very like the original relativistic 
case of / and Section 3 for which the results of Halperin are directly applicable. That is, W a b is diagonal, 



W ab (\r-r'\;t,t) = 6 ab W(\v-r'\;t,t), 



(54) 



whence W(\r — r' | ; t) = (</> a (r, t)(f> a (r' , t)) (no summation), the thermal Wightman function for either 0i or 02- 

For all the utility of bringing the effective potential into the action from the boundary conditions, W(r; t) is still 
most easily built from the modes U^ k , satisfying the equations of motion 



d 2 l2 



m 2 (t) 



0. 



(55) 



for 77i 2 (t) above, subject to the initial condition of a thermal distribution with chemical potential /x See the Appendix 
for greater detail. 

But for the chemical potential, this situation of inverted harmonic oscillators was studied many years ago by Guth 
and Pi [ffil and Weinberg and Wu pl[ | . In the context of domain formation, we refer to the recent work of Boyanovsky 
et al., [|22|| . For the case in hand, if we make a separation into the unstable long wavelength modes, for which |k| < M , 
and the short wavelength modes |k| > M, then W(r; t) is the real quantity 



W(r;t) = / fke lk - x C{k;fi) 

J\k\<M 

+ [ fke lkx C(k;fi) 

J\k\>M 



1 + A(k)(cosh(2Q(k)t) - 1) 
1 + a(k)(cos(2w(k)t) - 1) 



(56) 



with r = Ixl and 



Q 2 (k) = M 2 - |k| 2 
w 2 (k) = -M 2 + |k| 



A(k) 
a{k) 

All the /i-dependence is contained in the factor 

1 



2\ n 2 (k) 

if um 

2\ w 2 (k) 



2w(fc) 
where 

u 2 {k) = |k| 2 + m2 

In the zero chemical potential limit fi = mo + /iq — 0, C(k, /io) takes the familiar form 



(57) 



(58) 



(59) 
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CO, Mo) = -^-rcoth(/3 co(k)/2) (60) 

On the other hand, in the non-relativistic limit e(k) = k 2 /2m ~ lo — m and fio <C m C(fc,/io) is the equally familiar 
Bose distribution 

° {k > Mo) ~ 2^fc) ( l-e-^W-Mo) ) (61) 

Further, in the high temperature relativistic limit To ^> mo, C(k, /iq) of ( |60| ) simplifies as 

^co^(*)/2) * (62) 

Equally, in the non-relativistic limit when e(k) — <C Tq, 

C(k,fi )~ T ° — | (63) 
e(fc) + \fio\ 



from (p31h. This reproduces (K3 



C(fc, Mo )~_^ (64) 



k 2 + 



m 



up to an irrelavant coefficient of proportionality, on using the definition of mo given earler in (p2[). That is, in these 
regimes the relativistic and non-relativistic W(r, t) are identical, once we take the identifications of (p2]) into account. 

We observe that, if we were to use W(r;i) of ( |56| ) as it stands, then both W(0;i) and W"(0;i) necessarily suffer 
from ultraviolet divergences. However, the string thickness at the end of the quench will be 0{M~ 1 ). It is the zeroes 
coarse-grained to this scale that will provide the subsequent network. Thus, if we do not probe the field zeroes within 
a string we need consider only the first term in (^6|). There is another point. Even before the quench begins there is 
a high density of line zeroes coarse-grained to this same scale 0(M~ l ) in the initial equilibrium phase. From (|l^)) 
their density is ^] n(t) — 0(M 2 ) if M > m, which we assume. However, these modes are entirely transient due to 
the uncertainty principle. If we were to calculate the correlations of p^(x) at different times t and t' , we would find 
rapidly oscillating behaviuor with period At — 0(m _1 ). On the other hand, a calculation of the density correlations 
at different times from the unstable modes in (|56j ) does not give oscillatory, but damped, behaviour. It is the residue 
of the strings produced by the unstable modes that survives to produce the network, and the transient strings can be 
ignored. Henceforth, we retain only the first term 



W(r;t) = / 4 3 ke lkx C(k,vo) 

J\k\<M 



1 + A(k)(cosh(2Q(k)t) - 1) 



(65) 



of (p36|). In the two critical regimes where ( |62[ ) and (Q) are valid, W(r;t) is again the same for both cases. 

Even though the approximation is only valid for small times, there is a regime Mt > 1, for small couplings, in 
which t is large enough for cosh(2M£) w | exp(2Mi) and yet Mt is still smaller than the time (9(ln 1/A) at which the 
fluctuations begin to sample the ground-state manifold. In this regime 

W{r;t)~ f fkC{k 1 n)e ik - x A(k)e 2n(k)t (66) 
In these circumstances the integral at time t is dominated by a peak in the integrand fc 2 e 2n ( fe )* at k around k e , where 

tk *= M { 1+0 {w)} ^ 

and we have assumed M and mo to be comparable. The effect of changing (3q is only visible in the 0(1 /Mt) term. In 
fact, we are being unnecessarily restrictive in wanting to preserve the idnetical behaviour of (^2|) and ([34]). From (|l4| ) 



5 This was essentially the basis for Halperin's results on string densities in 112 
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onwards it follows that the overall scale of W is immaterial to the vortex density. All that is required for identical 
leading behaviour in relativistic and non-relativistic regimes is that C(fc,/Ko) varies slowly in the vicinity of the peak 
of the integrand at tk 2 ~ M. [A(k) is already slowly varying]. This is the case when, allowing for coefficients O(l), 



1< fir = —T < tM < In - (68) 
IMo| m 2 W 

where, in the same spirit, we have taken |//q| < The upper bound on tM is a reminder that interactions are 
always present, and the Gaussian approximation must fail as soon as the field fluctuations have extended to the true 
ground-states at the minima of the potential. The lower bound is necessary for the integrand to be peaked strongly so 
that the saddle-point approximation is valid. As long as ([38]) is basically correct, any difference between the relativistic 
and non-relativistic regimes will be non-leading. 

Assuming these limits, we recover what would have been our first naive guess for a correlation function based on 
the growth of the unstable modes, 

W{r-t)~ f fke lW - T e 2n{k) \ (69) 

J\k\<M 

We understand the dominance of wavevectors at k 2 = M/t in the integrand as indicating the formation of domains 
of mean size 



^{t) = 0(y/t/M) (70) 



once Mt > 1. Specifically, we take Q = 2y/t/M. Once Mt > 1 then £(f) > M _1 , where M~ x characterises 
the cold vortex radius. In the weak coupling approximation individual domains become large enough to accomodate 
many vortices before the approximation breaks down. There is no difficulty with causality since domains increase in 
size as £ = j= < 1. On neglecting terms exponentially small in Mt, W(r;t) of ( p9| ) can be further rewritten as 



W{r; t) ~ e 2Mt / dk sinc(fcr) k 2 e~ tk /M (71) 



,.2 



W(0;t)ex P [--^) (72) 



where 



e 2Mt 



WM * C WW' (73) 

for some C. The exponential growth of W(0;t) in t reflects the way the field amplitudes fall off the hill centred at 
$ = 0. With the peaking in wavelength I = fc _1 understood as indicating the appearence of domains of characteristic 
linear dimension the Gaussian in r is a reflection of the rms variation A£ in domain size £. This variation is 
large. If we isolate the Gaussian saddle-point in ( |69|) as 

W(r; t) ~ e 2Mt / dk sinc(fcr) k ^(k-k c ) 2 /2(Akf (74) 



then 



To calculate the number density of vortices at early times we insert the expression (|7lj) for W into the equations 
derived earlier, to find 

»(*) = ( 76 ) 



6 This somewhat arbitrary choice differs from that of |22] and I by a factor of y2. 
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for an 0(2) theory with strings in three dimensions. We note that the dependence on time t of both the density and 
density correlations is only through the correlation length We have a scaling solution in which, as the domains 
of coherent field form and expand, the interstring distance grows accordingly. Since the only way the defect density 
can decrease without the background space-time expanding is by defect-antidefect annihilation, we deduce that the 
coalescence of domains proceeds by the annihilation of small loops of string. However, because the density of vortices 
only depends on £(t) in this early stage, the fraction of string in 'infinite' string remains constant. Thus, at the same 
time as small loops disappear, other loops must rearrange themselves so that the length of 'infinite string decreases 
accordingly. Finally, there is roughly one string zero per coherence area, a long held belief for whatever mechanism. 

There is one final concern. A necessary condition for this rolling down of the field to b e valid is that the initial field 
distribution at t < should not overhang the point of inflection at (|^|) = -^/2M 2 /3A = 0(T C ). That is, the initial 
field fluctuations about <\> a = should be small enough that there is no significant probability that the field is already 
in the true vacuum. To check this for the relativistic regime, we anticipate that, when a domain structure forms, 
the smallest domains that can be identified will be of the size of the cold vortex radius, the Higgs field Compton 
wavelength, £o = 0(M^ 1 ). For temperature T > T c (but not too close), the rms field fluctuations on this scale are 
p7{ Aip — 0((TM)i). The condition that (A0) 2 < T 2 is guaranteed when T c is much larger than M, as happens 
for small coupling. In fact, since there are small prefactors, the coupling does not have to be very small for this to 
happen. Details are given in /. For the non-relativistic case there are no reasons, a priori, why the situation should 
be different. 



VI. VORTEX DENSITY CORRELATIONS 



In addition to the gross vortex density ( |76| ) we can calculate the density-density correlation functions C'ij(r;t) of 
(|l6|), identical for both the relativistic plasma and the non-relativistic medium. These were not considered in I, 

Yet again, as in the case of the density n(t), the i-dependence of Cy only occurs implicitly through The simple 
analytic form of (fnl) enables us to calculate A and B, up to exponentially small terms in Mt, as 



A(r;t) = 
B(r;t) 



7T 2 £ 4 (<) (l- e -2rV« 2 (t))2 
2 



(1 



em 



<o, 



5 -2r 2 /£ 2 (i) 



> o, 



tt 2 £ 4 W (l-e- 21 - 2 /^*)) 
whence the diagonal elements Ca (no summation) are all negative. Specifically, 

,,2 



2 / 2r 2 \ e - 2j - /« (*) 



(77) 



(78) 



That is, we have anticorrelation of densities for parallel directions and positive B for orthogonal ones. This is just 
the situation discussed in Section 2. 

It is useful to expand (|77| ) and ((7^) for r < £. On normalising by a factor of n 2 , a and b of (HI) are given by 



(79) 



a{r;t) =-1 + — < 



6(r;t) 



+ 0(1) > 0. 



For r > £ there is exponential falloff but, from ( [791) we see that, in units of (7r£ 2 ) -2 , the anticorrelation is large. The 
discussion of Section 2 is directly relevant. Since strings with a long persistence length would imply positive parallel 
correlations, and weak orthogonal correlations, we can interpret these anticorrelations as a reflection of an increased 
string bendiness. Although it is difficult to be precise, this suggests a significant amount of string in small loops. This 
is an important issue since, as we noted, early universe cosmology requires infinite string if string is to be the source 
of large-scale structure. In practice, some is enough. 

There is another, indirect, way in which we can see the tendency for more string to be in small loops than we might 
have thought. In numerical simulations of string networks based in throwing down field phases at random, the rule 
of thumb for a regular domain structure in field phase is that a significant fraction of string, if not most, is infinite 
string pLpl. 
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However, as we saw earlier in fl7q), we do not have a regular domain structure in our model but have domains with a 
large variance in their size A£/£ = i. Unfortunately, the domain structure that we have here is not yet appropriate for 
a direct comparison since, as well as the domain size, the field magnitude has a variance about (\<fi\) — 0(M e Mt ). The 
work of Guth and Pi shows that this can be parametrised by an effective dispersion in t in of At = 0(M~ 1 ). 
Nonetheless, consider an idealised case in which domain growth stops instantaneously because of back-reaction at 
some time tf. The distribution of strings will then be as above for £ = £(t/), while the field adjusts to the vacuum 
manifold, without changing phase, in each domain. We would expect some string- antistring annihilation to continue 
while this adjustment occurs, so that the n(t) calculated previously is an overestimate of the string density at the 
end of the transition. However, the domains will still be of varying size with variance plausibly given approximately 
by (|75|). Although it has not been introduced along the lines above, the inclusion of domain variance in numerical 
simulation of string networks shows p9| ] that, the greater the variance, the more string is in small loops. Beyond 
observing that there seems to be some infinite string, we will say no more. 

All our results are for quenches that go from significantly above the transition to significantly below it. These give 
the smallest value of A£/£ possible and, plausibly, the best chance of producing infinite string. Suppose the initial 
state characterised by So of (^6|) is very close to the transition. For relativistic theories this means starting from a 
temperature only just above the critical temperature. For a non-relativistic theory it corresponds to beginning from 
a density very close to the critical density. In either case it corresponds to taking TOq suffficicntly small that the 
inequalities (pq) cannot be satisfied. Some caution is required, since fluctuations are large near the transition and 
the semiclassical action So of(Mfl) is unreliable. Despite that, let us take it seriously and assume that it is possible to 
neglect m\ in comparision to fcj within the time interval 1 < Alt < ln(l/A). 

The overhang calculation of / shows that there is some leeway. Then, with 



T _ T 

+ TTln k 2 



(80) 



in (64) and (|62|) no longer slowly varying, instead of (69), we have 



W{r-t)~ / ^ e *k.x e 2n (fc) t (81) 

'|k|<M k 



dk sinc(fcr) e 2 ° (fc)t (82) 
up to irrelevant factors. 

There is now no peaking of the integrand in ( |82] ) at any preferential wave number, and hence no preferred length 
scale that can be identified as a typical domain size. This lack of domain structure, and the inapplicability of the 
Kibble mechanism, is because of the enhanced long-wavelength fluctuations that come from being too close to the 
transition initially. We stress that this does not mean that there are no vortices, or even a small vortex density. 
W(r; t) of ( |S2"| ) can be integrated as 

W(r;t) = W(0 ] t) 1 F 1 (±l,--£-) (83) 



.2' 2' 

for £(i) as before. A little algebra shows that the vortex density n{t) is 

one-third of its previous value in (|7(]). For distances r > £, W(r; t) now has power behaviour 

W(r;t) = o(^j (85) 

instead of the Gaussian behaviour of (f72|). As before, there is anticorrelation with A < 0, B > but, at large distances, 
the density-density correlation functions are very different, only showing power-behaviour in their fall-off. However, 
for r ~ £ their behaviour is very like that of (|7^). Nonetheless, if we interpret ( jB2|) as the limiting case of domain 
structure with maximum variance in size (compatible with causality), the same numerical data quoted earlier p9[ , 
with due caveats, suggest even more string in small loops. 



Finally, with superfluid films in mind, we can repeat the analysis in two space dimensions. Instead of (72) W(r;t) 
is given, when Mt > 1, by 
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W(r- 1) ~ e 2Mt / dk sinc(fcr) fc e - tfe2/M (86) 



The integrand is strongly peaked at k 2 = M/2t which, from our previous analysis, we interpret as the existence of 
domains, trapping 0(2) monopoles (vortex cross-sections) on their boundaries. The domains are of varying size. If 
we isolate the Gaussian saddle-point in (R5I) as 



W(r; t) ~ e 2Mt I dk sinc(fer) kce -(k-k c f /2(Akf (8?) 



then the rms variation A£ in domain size £ is larger than in three dimensions, as 

A£ _ Afc _ 1 

This is independent of the definition of £, but if we take £ 2 (t) = 4/fc 2 as before, then 

3 2r 2 

~e(t) 



W(r;t)=W(0;t) 1 F 1 [l,-,- 7 ^) (89) 



for r > £. This is intermediate between ( p5| ) and the Gaussian behaviour of (|72j). Yet again there is anticorrelation in 
the density-density correlation function (p(r)p(O)) — A(r;t), where A(r;t) is given in (|l7|). 



VII. CONCLUSIONS 



In this paper we have shown how global 0(2) vortices appear, at a quench from the ordered to disordered state, 
as a consequence of the growth of unstable Gaussian long wavelength fluctuations. Most importantly in the light 
of discussions about the extent to which vortex production in low-temperature many-body systems simulates vortex 
production in the early universe, our model supports the analogy. We have shown how, with our simple assumptions, 
vortex production is identical in both a relativistic high-temperature quench, in which the initial state is characterised 
by T > m, and in a non-relativistic density quench in which the initial state is described by T <C m (to freeze out 
antiparticles) with a chemical potential fi nr ~ T . All that is required is an appropriate translation of the parameters 
from the one case to the other. 

In our simple model of Gaussian fluctuations the resulting string configurations scale as a function of the correlation 
length £(t) = 0(t 3 ), at about one vortex/correlation area. This is compatible with the Kibble mechanism for vortex 
production on domain boundaries upon phase separation. For a weak coupling theory the domain cross-sections are 
significantly larger than a vortex cross-section at the largest times for which the approximations are valid. However, 
there is a large variance in their size, with A£/£ = i. Because of this there is more string in small loops than we 
might have anticipated. Nonetheless, we expect some infinite string. 

We stress that our model can, at best, describe weak coupling systems for the short times while the domains 
are growing before the defects freeze out. This is unsatisfactory for most early universe applications and for low- 
temperature many-body systems. However, we know in principle [p2[ how to include back-reaction (still within the 
context of a Gaussian approximation) to slow down domain growth as the field fluctuations spread to the ground-state 
manifold. The identity of the weak-coupling results of the two regimes should survive to this case also, although it 
will probably lead to different conclusions from those above. Further, there is no difficulty in principle of embedding 
these results in a FRW metric, along the lines of pPf . 
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APPENDIX 

We choose to work with 'physical' fields, that is fields which evolve according to the original Hamiltonian, Hq, which 
has not had the chemical potential absorbed into it. The chemical potential then appears in the boundary conditions 
and so is placed on the same footing as the temperature. For this reason we work with the eigenstates of Q, 4>(x) and 
${x), rather than the Cartesian components (f>i(x) and <f>2{x). The conversion to the Cartesian thermal Wightman 
functions W a b(r,t) of the text is trivial. 

For the complex scalar field we are interested in, the Hamiltonian has the usual form as given in (|3^) in which, for 
the reasons given in text, A = 0. That is, we are interested in 'free' thermal Wightman functions, albeit with changing 
mass. The Wightman functions in question are defined to be 

G>(x,x') = Tr{e-^ H °-^(l){x)^(x')} (Al) 
G<(x,x') = Tr{e-^ H °-^^(x')(p(x)}. (A2) 

Other types of propagator are defined in a similar way, retarded, advanced, time-ordered etc. but will not be needed 
here and in any case they can be built from G >( - < ^ relatively easily. 
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The free propagators satisfy differential equations which look like the usual Klein-Gordon equation, the only differ- 
ence in our case is that the mass term is time dependent - a simple step function for the quench at time t — that we 
have considered in the text. Otherwise our dicsussion is general. Adopting it to our specific model is straightforward. 

We start by solving the homogeneous equation 

{d 2 t +uo 2 {t))U ± (t) = (A3) 

where 

u(t) = 6{-t)wi +6{t)w 2 . (A4) 

As it stands both u}\ and u>i are assumed to be real. If we have an unstable mode at late times then we just replace 
uo 2 = ±ifi. The boundary conditions are 

lim U ± {t)=&c£,{±iu 1 t} (A5) 

t — ► — CO 

U ± (t = + ) = U ± (t = 0-) d t U ± {t = Q + ) = d t U ± {t = Q-). (A6) 

This gives 

f7±(i) = 8{-t)e ±Wlt + 6{t)[A±e ±lUlt + B ± e^ tuJlt } (A7) 

Now we can construct the thermal Wightman functions. As they satisfy 

{d 2 +u 2 {t))G >{<) {t 1 t') = (A9) 

(Of, +Lu 2 (t'))G >{<) (t,t / ) = (A10) 

we try 

G >{<) (t,t / )^a >(< '>U-(t')U + (t)+d > ^U + (t')U-(t) (All) 

where terms such as U~(t')U~(t) and U + (t / )U + (t) are eliminated by demanding time translation at early times. 

We start by noting that it is common, when considering the thermal Wightman functions of real fields in thermal 
equilibrium, to use their definitions to deduce a trivial identity relating the two thermal Wightman functions at equal 
times. In the case of the complex fields we see that we have 

[G >(<) (M')]* = G >(<) (t'*,r)G >(<) (M) - G >(<) (i*,0- (A12) 

We have been very careful to ensure that the definition of ${t) is such that it evolves like <fr(t), namely ${t) — 
exp{— iHot}(jJ (0) exp{ii?oO (where we have dropped spatial labels). Problems occur with complex time shifts, such 
as are encountered in thermal field theory, as roughly speaking we require (</>(£))* = (jr (t*). The U + and U~ functions 
below are related in a similar way. 

The ETCR (equal time commutation relations) can be used to provide boundary conditions, 

ld>(t),W)} = l<PHt),nHt)]=i- (Ai3) 

We have U(t) = d t (f>^(t) and 11+ (t) = d t <p{t) since, for the physical fields, there is no chemical potential in the 
expressions for IT and 11+ . Looking at these equations for any one negative time gives 

a< -d< -a> +d> = — (A14) 

and we find the same relation if the ETCR is imposed at any positive time. This confirms that the ETCR are indeed 
maintained at all times if they are true at any one time. 
There is a second ETCR which we can use; 

W)A\t)]=Q (A15) 

This requires 
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G>(t,t) = G < (t,t) (A16) 

which gives 

a > +d > =a < +d < . (A17) 

Once this is satisfied it ensures that this ETCR also holds at all times. 

The temperature and chemical potential appear only in the last boundary condition. This is the well known KMS 
condition |3l| . Here we enforce this at any negative time as then the system is in equilibrium. Thus 

G>(t,t') = e-^G<(t + i/3,t') (A18) 

This holds for all t, t' < 0, but not for later times, so we find 

a> = e- (3 ^ 1+ ^a < (A19) 
d> = e l3 ^ 1 -^d < (A20) 

Putting all this together gives 

(A21) 
(A22) 

d< = TT 7>T 7 (A23) 

d> = 7T h v (A24) 

For our purposes we just need the equa l tim e propagator which can be obtained from either of the thermal Wightman 
functions because of the second ETCR (A16). From the above we find that 

i e~^ Wl — e^ 1 

G > (t,t) = G < (t,t) = s 3 a --.U-(t)U + (t) (A25) 

from which our results of section 5 follow. 

We might ask whether we need worry about off-diagonal terms in the Cartesian propagator or, equivalently, 
(cf>(t)<fi(t')) and (0t (t)$ (t')) . In some approaches using different definitions of the fields, e.g. p5| , these are not 
zero. Trying the same ansatz as before, we see that time-translation invariance and the KMS condition, applied at 
negative times, ensure that these are zero for all times in our case. This is a big advantage of our definitions. 
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